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.expected,  whereas  for  lower  incidence  angles  the  base  cavity  and  sometimes  double 
cavities  appeared.  The  trends  of  the  force  coefficients  for  these  different  cavity- 
flow  patterns  were,  of  course,  found  to  be  totally  different.  This  fact,  which  was 

never  considered  in  the  design  procedure  nor  in  the  off -design  prediction  theory, 

might  have  caused  an  erroneous  prediction  for  the  thrust  coefficient  and  efficiency  j 
of  this  supercavitating  propeller. 

In  order  to  make  comparisons  with  the  experimental  data,  three  new  nonlinear  cavity 
flow  theories  have  been  developed,  one  for  the  supercavitating  (S/C)  regime  and  two  j 
for  the  partially  cavitating  (P/C)  regime.  The  results  of  these  theories  compared 
favorably  with  the  experimental  data,  but  the  accuracy  degraded  as  the  cavity  length  , 
became  close  to  the  chord  length  both  in  S/C  and  P/C  ranges.  The  success  of  the 

present  experiments  has  allowed  us  to  design  dummy  blades  for  the  cascade  experiments 

to  be  conducted  in  the  following  phase  of  the  present  work  for  FY  1981. 
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ABSTRACT 


I 


Experiments  were  conducted  at  the  High  Speed  Water  Tunnel  of 
the  California  Institute  of  Technology  for  measuring  the  flow 
characteristics  of  a  single  foil  over  a  full  range  of  cavita¬ 
tion  numbers.  The  blade  profile  shape  used  is  that  of  a 
modified  Tulin  two-term  camber  with  a  blunt  trailing  edge. 

This  profile  shape  was  taken  after  the  cross-section  profile 
of  a  supercavitating  propeller  designed  for  an  actual  high 
speed  hydrofoil  boat.  It  was  discovered  in  the  experiments 
that  +-he  flow  pattern  totally  changed  at  about  an  incidence 
angle  of  2  degrees.  For  higher  incidence  angles  than  2  degrees, 
the  cavity  sprang  from  the  leading  edge  as  was  expected, 
whereas  for  lower  incidence  angles  the  base  cavity  and  some¬ 
times  double  cavities  appeared.  The  trends  of  the  force  coef¬ 
ficients  for  these  different  cavity-flow  patterns  were,  of 
course,  found  to  be  totally  different.  This  fact,  which  was 
never  considered  in  the  design  procedure  nor  in  the  off-design 
prediction  theory,  might  have  caused  an  erroneous  prediction 
for  the  thrust  coefficient  and  efficiency  of  this  supercavitating 
propeller . 

In  order  to  make  comparisons  with  the  experimental  data,  three 
new  nonlinear  cavity  flow  theories  have  been  developed,  one 
for  the  supercavitating  (S/C)  regime  and  two  for  the  partially 
cavitating  (P/C)  regime.  The  results  of  these  theories  com¬ 
pared  favorably  with  the  experimental  data,  but  the  accuracy 
degraded  as  the  cavity  length  became  close  to  the  chord  length 
both  in  S/C  and  P/C  ranges.  The  success  of  the  present  experi¬ 
ments  has  allowed  us  to  design  dummy  blades  for  the  cascade 
experiments  to  be  conducted  in  the  following  phase  of  the 
present  work  for  FY  1981. 
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1.0  INTRODUCTION 


Recently  a  theory  for  predicting  the  off-design  performance 
of  cavitating  propellers  was  developed  by  Furuya  (1976,  1978) 
and  Furuya  and  Maekawa  (1979)  based  on  the  supercavitating 
and  partially  cavitating  cascade  theories.  The  results  of  the 
theory  compared  favorably  with  experimental  data  except  for 
the  large  J-range  at  which  the  partially  cavitating  condition 
prevails  over  the  propeller  blades.  A  question  did  arise  as 
to  the  accuracy  of  the  supercavitating  and  partially  cavitating 
cascade  theories  of  Furuya  (1975)  and  Furuya  and  Maekawa  (1979) 
which  were  used  to  form  the  sectional  loading  data  in  the  above 
propeller  theory.  Only  one  experimental  work  for  the  cavitating 
cascade  exists  to  date.  The  work  was  conducted  by  Wade  and 
Acosta  (1969)  who  used  a  plano-convex  cascade  blade,  which  is 
quite  different  from  practical  blade  profiles  used  for  high 
speed  propellers. 

The  objective  of  the  present  study  is,  therefore,  to  conduct 
water  tunnel  experiments  for  cavitating  cascades  having  prac¬ 
tical  blade  profile  shapes  and  to  compare  the  results  with  the 
existing  theories.  The  program  was  divided  into  two  phases, 
Phase  A  and  Phase  B.  We  just  completed  Phase  A  for  FY  1980, 
the  results  of  which  are  reported  herein.  This  phase  is  what 
is  called  a  preparatory  phase  to  the  cascade  experiment  of 
Phase  B  for  FY  1981.  Before  constructing  an  experimental 
set-up  for  cascade,  it  was  imperative  for  us  to  conduct  a 
single  foil  experiment  to  ensure  the  appropriateness  of  the 
blade  profile  and  dimensions.  We  chose  a  cross-sectional 
profile  shape  of  Hydronautics  7607.02  supercavitating  pro¬ 
peller  at  50%  radial  station.  This  cross  section  profile  is 
made  of  a  modified  Tulin  two-term  camber  with  a  blur.c  trailing 
edge  (see  the  report  of  Furuya  (1978)  for  the  modifications 
in  detail) . 


We  fabricated  a  foil  based  on  the  above  off-set  data  and  con¬ 
ducted  single-foil  experiments  at  the  High  Speed  Water  Tunnel 
at  the  California  Institute  of  Technology.  The  data  were  taken 
over  a  wide  range  of  cavitation  numbers  covering  from  the 
fully-wetted,  partially  cavitating  and  supercavitating  regimes. 
Although  the  purpose  of  the  program  was  to  compare  theories 
with  experiments  for  cascade  flows,  it  was  considered  that  a 
similar  comparison  for  the  single  foil  data  just  obtained  would 
be  interesting.  Despite  many  years  of  studies  for  cavity  flows 
to  date,  it  was  a  surprise  to  learn  that  there  existed  no 
appropriate  nonlinear  theory  for  determining  the  characteristics 
of  cavitating  foils  having  practical  foil  profile  shapes  in 
the  unbounded  flow  media.  Wu  and  Wang  (1964)  applied  an  open- 
wake  model  for  supercavitating  foils,  but  the  results  for  the 
drag  force  prediction  were  not  quite  satisfactory  due  to  the 
choice  of  the  flow  model.  Larock  and  Street  (1965,  1968)  used  a 
more  reliable  wake  model,  i.e.,  Tulin's  single  vortex  wake  model 
to  the  S/C  foil,  but  their  solution  method  was  applied  either  to 
flat-plate  foils  or  an  inverse  type  of  camber  specification. 

For  the  partially  cavitating  regime,  not  even  a  single  nonlinear 
theory  existed. 

It  was  for  this  reason  that  we  decided  to  devote  our  effort 
to  the  development  of  nonlinear  partially  cavitating  and  super¬ 
cavitating  theories  which  could  be  applied  to  any  practical 
foil  profile  cases.  With  these  theories  developed,  it  was 
believed  that  we  could  fill  the  gap  existing  in  the  cavity 
flow  theory  regime.  For  the  supercavitating  flow,  the  theory 
of  Larock  and  Street  (1965)  who  applied  the  single  vortex 
model  was  extended  for  general  blade  profile  cases.  On  the 
other  hand,  we  developed  two  partially  cavitating  flow  theories, 
one  using  the  double  wake  model  and  the  other  using  the  open 
wake  model. 

It  must  be  pointed  out  here  that  the  results  of  the  present 
work  will  form  the  basic  foundation  for  conducting  the  cas¬ 
cade  work  in  Phase  B.  Since  it  was  found  in  the  present 


t 


experiment  that  the  blade  profile  shape  and  its  dimensions  were 
appropriate,  dummy  blades  to  form  a  cascade  configuration  were 
already  designed.  Four  dummy  blades  will  be  fabricated  in  the 
following  phase.  Furthermore,  the  foil  used  for  the  present 
experiment  will  be  used  as  one  of  the  cascade  blades,  i.e., 
the  center  foil  for  the  force  measurements. 

In  what  follows,  we  will  describe  the  experimental  work.,  theory 
development,  and  data  comparison  between  the  theories  and 
experiment . 


2.0 


EXPERIMENT 


2.1  MODEL 

The  final  goal  of  the  program  is  to  obtain  the  experimental 
data  for  cavitating  cascades  having  practical  blade  profiles 
and  to  make  comparisons  with  the  existing  nonlinear  theories. 
As  has  been  mentioned  in  the  preceding  section,  due  to  the 
limited  time  and  funding  we  decided  to  devote  the  first  phase 
of  the  program  mainly  to  a  single-foil  study  which  will  form 
a  basis  for  the  study  of  the  cascade  configuration  to  be  con¬ 
ducted  in  the  following  phase. 

The  difference  between  the  present  work  and  numerous  previous 
works  lies  in  the  fact  that  the  present  one  investigated  the 
hydrodynamic  characteristics  for  foils  having  practical  blade 
profile  shapes.  By  "practical"  we  mean  that  the  blade  profile 
must  be  different  from  flat-plate  or  circular  arc  profile 
shapes.  For  achieving  the  purpose,  the  cross-section  profile 
of  Hydronautics  7607.02  S/C  propeller  at  50%  radial  station 
(r  =  .5)  was  chosen.  The  upper  and  lower  blade  off-set  data 
are  given  in  Table  2.1  (see  also  Bohn  and  Altman  (1976)).  The 
reasons  for  choosing  the  cross-section  profile  at  r  =  .5  were 
based  on  the  solidity  and  stagger  angle  of  the  propeller  which 
would  be  readily  implemented  in  the  High  Speed  Water  Tunnel 
(HSWT)  at  the  California  Institute  of  Technology.  In  order 
to  carry  out  the  cascade  experiment,  at  least  five  blades  will 
be  desired.  By  choosing  the  chord  length  of  the  blade  to  be 

3.2  inches,  the  high  solidity  (sol)  such  as  sol  =  1  will 
possibly  be  achieved.  Since  the  two-dimensional  test  section 
available  at  HSWT  was  made  to  have  the  stagger  angle  of  45°, 
that  of  the  Hydronautics  7607.02  S/C  propeller  at  r  =  .5  which 
is  48.9°  (see  Furuya  (1978)  for  more  details)  matches  well 
with  this  number.  As  far  as  the  practicality  of  the  blade 
profile  is  concerned,  the  profile  chosen  here  is  essentially 

a  Tulin's  two-term  camber  but  with  the  leading  edge  (x  <  0.2) 


and  trailing  edge  (x  >  0.8)  slightly  modified.  The  report  by 
Furuya  (1978)  describes  these  modifications  in  detail. 

The  blade  set-up  to  the  tunnel  wall  is  shown  in  Figure  2.1. 

Two  pressure  taps  were  placed  on  the  model.  One  was  placed 
at  20%  chord  from  the  leading  edge  on  the  suction  side  to 
measure  the  pressure  inside  the  cavity  and  the  other  at  the 
base  of  the  foil  to  measure  the  base  pressure  inside  the 
wake . 

Grooves  were  milled  in  the  upper  surface,  and  1/16  inch  outer 
diameter  brass  tubes  were  laid  in  these  grooves.  Transparent 
epoxy  resin  was  molded  into  the  grooves,  completely  covering 
the  brass  tubes,  and  then  ground- smoothed  with  the  curved 
plane  of  the  upper  surface  of  the  model.  Holes  of  0.020  inch 
diameter  were  drilled  at  the  specified  pressure  tap  locations. 

The  high  lifting  force  expected  at  the  partially  cavitating 
and  the  fully  wetted  condition  required  the  use  of  the  heat 
treated  stainless  steel  17-4PH.  The  model  was  fabricated  and 
measured  to  an  accuracy  of  +  0.002  inches.  The  leading  edge 
is  specified  to  be  sharp,  however,  sharpness  in  the  most 
physically  possible  sense. 

The  cascade  experiment  to  be  conducted  in  the  following  phase 
simply  requires  dummy  blades  which  will  be  mounted  to  the 
other  side  of  the  tunnel  test  section.  These  dummy  blades 
have  already  been  designed  (see  Figure  2.2)  and  are  ready  for 
fabrication.  Four  dummy  blades  with  the  center  blade  (just 
used  for  the  present  experiment)  will  form  the  cascade  con¬ 
figuration  with  a  solidity  of  unity.  However,  with  only  three 
blades  used,  i.e.,  two  dummy  blades  plus  the  center  blade,  we 
can  construct  a  cascade  having  a  solidity  of  0.5. 


2.2 


EXPERIMENTAL  SET-UP  AND  APPARATUS 


The  experiments  were  conducted  in  the  High  Speed  Water  Tunnel 
(HSWT)  of  the  Hydrodynamics  Laboratory  at  the  California  Institute 
of  Technology.  The  tunnel  is  a  closed  circuit  water  tunnel 
having  an  available  pump  power  of  approximately  350  hp  to 
drive  water  through  the  working  section  at  various  controlled 
pressures.  The  speed  of  the  water  flow  in  the  working  section 
of  the  tunnel  can  be  varied  from  a  very  low  value  of  a  foot 
per  second  up  to  the  maximum  of  80  ft/sec.  More  basic  features 
of  the  water  tunnel  are  discussed  at  length  in  Knapp  et.  al . 

(1948)  and  Ward  (1976)  . 

For  this  experiment,  we  used  the  two-dimensional  working  sec¬ 
tion,  the  dimensions  of  which  are  30  inches  high,  50  inches 
long,  and  6  inches  wide.  In  our  experiment,  the  two- 
dimensional  working  section  was  installed  by  rotating  90°, 
thus  horizontally,  to  avoid  a  large  static  pressure  gradient, 
a  particular  concern  for  the  upcoming  cascade  experiment. 

The  foil  was  mounted  to  a  base  fixture,  fitted  with  a  balance 
port  fairing  disk.  The  force  balance  available  for  the  HSWT 
has  a  low  and  a  high  force  range  capacity.  The  experiment 
used  the  low  range  load  cells  having  a  range  of  +  200  lbs. 
lift,  +  80  lbs.  of  drag  and  +  300  in-lb.  of  moment.  The 
accuracy  of  each  cell  is  within  +  0.25%  of  the  rated  load. 

This  accuracy  is  inclusive  of  all  hysteresis,  non-linearity, 
and  repeatability. 

The  force  balance  data  as  well  as  the  water  tunnel  flow 
velocity  and  pressure  were  stored  in  the  buffer  of  the 
electronic  data  acquiring  system  at  the  HSWT  and  punched  on 
regular  IBM  cards.  These  cards  were  conveniently  used  as 
input  data  for  data  reduction.  The  data  acquisition  system 
also  has  an  integration  capability  for  the  force  data  over  any 
time  duration  so  that  the  time  averaged  static  forces  can  be 
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recorded.  The  force  balance  accuracy  mentioned  before  included 
the  stability  of  all  the  electronics  used  in  this  system. 

The  tip  clearance,  which  is  the  space  between  the  tip  of  the 
foil  and  the  opposite  wall,  was  adjusted  to  approximately  0.07 
inches  to  avoid  wall  contact  when  the  tunnel  pressure  is 
reduced  to  the  minimum  value. 

2 . 3  PROCEDURE 

The  model  was  set  at  angles  of  attack  a  =  0,  2,  4,  6,  8,  and 
10°  at  zero  water  velocity.  Then  the  water  velocity  was  in¬ 
creased  to  =  20  ft/sec  or  30  ft/sec,  depending  upon  the 
angle  of  incidence.  Originally  all  the  runs  were  scheduled 
with  30  ft/sec,  however,  the  large  force  increase  at  the 
transition  region  between  the  partial  cavitating  and  super- 
cavitating  conditions  forced  us  to  reduce  the  water  velocity 
down  to  20  ft/sec,  particularly  at  the  higher  incidence 
angles,  a  =  8°  and  10°. 

The  pressure  in  the  working  section  was  then  gradually  de¬ 
creased  from  the  atmospheric  pressure  until  the  cavity  started 
appearing  at  the  leading  edge  of  the  blade  and  finally  choked 
the  tunnel.  The  lift  and  drag  data  were  taken  along  with 
the  top  and  side  photographs.  All  the  force  measurements 
were  time  averaged  over  2  seconds.  The  top  photographs  were 
taken  with  4x5  inch  professional  film  centered  at  approxi¬ 
mately  3.7  chord  length  downstream  of  the  model.  The  side 
photographs  were  taken  with  35-mm  film  centered  approximately 
3.6  chord  length  downstream  of  the  model.  Both  photographs 
were  synchronized  by  the  strobe  light.  The  lengths  of 
cavities  were  measured  from  these  photographs. 

The  water  speed  was  determined  from  the  pressure  difference 
between  the  test  section  and  the  constant  cross-sectional  area 
at  the  upstream  of  the  contracting  nozzle.  The  pressure 
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difference  was  measured  by  a  mercury  manometer.  The  static 
pressure  of  the  tunnel  test  section  was  also  measured  at  the 
top  of  the  section  by  means  of  a  mercury  manometer. 

The  vapor  pressure  of  water  was  used  for  the  calculation  of 
the  approximate  cavitation  number,  but  was  corrected  for  pre¬ 
sentation  of  the  experimental  data.  More  details  will  be 
discussed  in  the  following  section. 

2.4  DATA  REDUCTION 

Various  corrections  were  made  on  the  cavitation  numbers  and 
the  results  of  the  force  measurements.  The  cavitation  number 
based  on  the  cavity  pressure  is  defined 


where  p^  is  the  static  pressure  at  the  center  of  the  test 
section  (3  inches  from  the  top  and  3  inches  from  the  bottom), 
t  is  the  density  of  water,  is  the  velocity  of  the  water  at 
the  test  section,  and  pc  is  the  cavity  pressure.  On  the  other 
hand,  the  cavitation  number  based  on  the  water  vapor  pressure 
is  defined 


where  pv  is  the  vapor  pressure  of  water  at  the  temperature 
during  the  experiment.  Initially,  we  had  planned  to  measure 
the  pressure  inside  the  cavity  by  using  one  of  the  pressure 
taps  mounted  on  the  foil.  However,  we  had  difficulty  in  doing 
so,  particularly  for  short  cavity  cases.  Also,  when  the  tunnel 
pressure  was  reduced  to  small  values,  there  existed  an  air 
leakage  problem,  which  caused  an  error  in  measurement.  We  had 
a  very  short  tunnel  time  assigned  and  there  was  no  time  for 
repair . 
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It  was,  therefore,  decided  to  make  corrections  on  a  's  to  pro¬ 
vide  all  the  data  in  terms  of  a's.  The  correlation  between 
a  and  was  presented  by  Kermeen  (1956)  and  Wade  and  Acosta 
(1965)  and  this  was  used  for  the  present  purpose.  Some  of 
the  measured  cavity  pressure  data  were  utilized  for  justifying 
this  correction  method.  Due  to  the  lack  of  correction  data 
in  the  region  a  >  1.5,  the  cavitation  numbers  are  left  uncor¬ 
rected.  It  must  be  mentioned  that  the  use  of  uncorrected 
cavitation  numbers  for  large  a's  will  not  greatly  affect  the 
force  data  presentation  since  the  force  is  fairly  flat  as  a 
function  of  a. 


The  measured  forces  were  normalized 

r  =  Lift 
L  ^pq^A 


c 


D 


-Ora? 

J5Pq12A 


where  A  is  the  plan  form  area  of  the  model.  First  it  was 
found  that  the  influence  of  the  pressure  change  of  the  tunnel 
during  the  experiment  on  the  force  balance  reading  was  con¬ 
firmed  to  be  negligibly  small.  Also,  based  on  the  previous 
literature  such  as  that  of  Barker  and  Ward  (1976) ,  very  slight 
corrections  were  made  for  the  measured  lifting  forces.  Thus, 
only  the  drag  forces  were  corrected.  The  corrections  included 
the  viscous  drag  on  the  model  and  the  fairing  disk.  Over  the 
tested  water  velocity  range,  the  boundary  layer  is  considered 
to  be  predominantly  laminar.  The  viscous  drag  on  the  model 
and  fairing  disk  was  estimated  by  the  Blasius  boundary  layer 
equation.  The  wetted  areas  for  this  correction  are  varied 
depending  on  the  fully  wetted,  partially  cavitating,  or  super- 
cavitating . 


The  wall  effects  for  the  cavity  flows  in  solid  wall  tunnels 
are  discussed  by  Whitney  and  Brennen  (1969)  and  Baker  (1972). 
Since  the  method  of  Baker  (1972)  indicated  such  a  correction 
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to  be  small  for  our  tunnel  set-up,  the  data  were  left  uncor¬ 
rected  for  the  wall  effect. 

All  other  corrections,  such  as  buoyancy,  wall  boundary  layer 
interference,  solid  blockage  and  wake  blockage  were  not 
included . 
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3.0  THEORY  DEVELOPMENT 

3.1  SUPERCAVITATING  FLOW  THEORY 

Figure  3.1  shows  a  flow  configuration  of  a  two-dimensional 

supercavitating  foil.  In  order  to  represent  the  cavity 

termination  and  highly  turbulent  flow  subsequent  to  it,  the 

single  spiral  vortex  model  of  Tulin  (1964)  is  employed.  On 

the  blade  the  cavity  is  assumed  to  spring  from  the  sharp 

leading  edge  B  and  trailing  edge  T,  both  fixed  but  arbitrary, 

extended  downstream.  The  pressure  inside  the  cavity  is  taken 

to  be  constant,  p  .  The  blade  reference  coordinate  is  in- 

c 

dined  at  an  angle  of  a  in  an  undisturbed  flow  so  that  the 
incoming  uniform  flow  is  always  fixed  to  be  parallel  to  the 
x-axis  of  the  coordinate  system.  There  exists  a  reason  for 
not  inclining  the  flow  angle  but  instead  tilting  the  foil. 

The  reason  will  be  described  during  the  course  of  the  theory 
development.  The  origin  of  the  coordinate  is  placed  at  the 
leading  edge  of  the  blade  with  x-axis  parallel  to  the  undis¬ 
turbed  flow.  The  wetted  portion  of  the  blade  is  defined  by 
y  =  f(x)  or  by  its  inclination  3  =  tan  1 ^ ^ j .  The  physical 
flow  field  is  mapped  onto  the  potential  plane  w  as  shown  in 
Figure  3.2.  The  complex  potential  of  the  stagnation  point  is 
chosen  to  be  (0,0)  in  this  transformation.  The  potential 
plane  is  then  transformed  onto  the  upper  half  of  a  new  plane 
;  =  '  +  in  of  Figure  3.3  by  a  mapping  function 


;  a\Z$~-"w  (3-1) 

or 

W  =  - 5-  (3.2) 

a  +  $  . 

The  cavity  end  point  is  mapped  to  infinity  and  the  infinity 
point  in  z  or  W  plane  is  now  mapped  onto  a  point  ;  =  ia. 
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A  hodograph  variable  w  is  introduced 


dW 

dz 


^ce 


-iw 


where  q  and  9  are  the  magnitude  and  direction  of  the  flow 
velocity.  Thus 


oj  =  9  +  ix 


where 


and 


where  the  uniform  flow  velocity  q1  is  taken  to  be  unity 
throughout  the  present  analysis. 


The  boundary  conditions  on  the  real-axis  £  are  now  expressed 
either  in  terms  of  9  or  x; 


(i)  x  =  0 

(ii)  9=8 

(iii)  9  =  it  +  8 


-oo  <  £  <  -  1  and  b  <  £  <  °° , 

-  1 <£  <0, 

0  <  £  <  b . 


This  is  a  typical  mixed-type  boundary  value  problem,  and  the 
solution  for  w  is  readily  written 


wU)  =  -y/(S  +  !)(?-  b) 


_ 26 _ 

i  yfU  +  £  ' )  (b  -  £  ' ) 


+ 


1 

2  xr  i 


b 


/ 


_ 2tt _ 

i^(l  +  £’)  (b  -  £’) 


(3.3) 
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(ii) 


i 


f  1  =  Re  <  oj  (ia)>  =  0 


f  2  =  Im  (  u)  (ia)\  -  In 


(3.4) 

(3.5) 


wetted  arc  length  =  S : 

On  the  wetted  part  of  the  streamline, 


dz  _  iS 
ds  e 


or 


ds  = 


-iS  e  dW 

6  5; 


(3.6) 


For  -1  <C  <  b,  oj(C)  can  be  written  as  follows: 


'D(')  = 


|i  g(?)  +  3(C)  /  -1<5  <  0 

1  i  g ( C)  +  BU)  +  *  ,  0  <5  <  b 


where 


9(5)  =  yj  (1  +  O  (b~5) 


V 

- 1  / _ iio _ 4^ 

T  /  yj(l  +  C  1  )(b  -  £  ~)  ’ 


_  f  _ 1  _  dS  ^ 

J  yj(l  +  C  ')(b  -  •  ')  ’  '  "  5 


-7-  +  P 


(3.7) 


and  j  indicates  a  Cauchy  principle  integral.  By 
integrating  Equation  (3.6)  we  obtain 


s  ( 


b 

-■/ 


sgm  *--■  d: 

Vl  +  o 


(3.8) 


i 

j 

i 
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where 
sg(5)  = 


i  1  C  >  0 

|  9 

-1  5  <0 


The  arc  length  condition  is  therefore  satisfied  by 
the  following  equation: 


f3  5  S  -  s(-l)  =  0 


(3.9) 


(iii)  closure  condition: 

For  zero  net  source  strength  within  the  contour  C 
enclosing  the  body-cavity  system  this  can  be  written 


'/■ 


Re  <f  t  w  (  S )  dW  )  =  0. 
C 


(3.10) 


And,  after  transformation  (see  p.  173  of  Larock  and 
Street  (1965)), 


Re 


j/-<‘>  S?«j- 


(3.11) 


where 

dW 

d? 


2a  »Lc 
U2  +  a2,2  ' 


Since  C  has  a  double  pole  at  £ 
dition  becomes 


=  ia,  the  closure  con- 


or 


f  A  =  Re 


(ia)>  =  Re 


H' (ia) 
H  (ia) 


oj(ia)  +H(ia)  G'(ia)}  =  0  ( 


where 

H  =  yjTi  +  1)U  ~  b) 


.12) 


T 
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and  G  is  the  remaining  part  of  the  right  hand  side  of 
Equation  (3.3) . 


A  set  of  four  nonlinear  equations  just  obtained,  i.e., 
(3.4),  (3.5),  (3.9)  and  (3.12),  will  determine  the  four 

unknown  solution  parameters. 


It  must  be  mentioned  here  that  when  one  applies  the 
transformation  formula  for  the  closure  condition  from 
Equation  (3.10)  to  Equation  (3.11),  it  will  become 


Re 


“iotl  jduU) 

I 


0  . 


Therefore,  if  there  exists  a  finite  incoming  flow  angle, 
the  closure  condition  takes  a  different  form  from 
Equation  (3.10).  This  is  the  reason  for  maintaining 
the  zero-incoming  flow  set-up  as  has  been  mentioned 
earlier . 


The  lift  and  drag  coefficients,  CT  and  C_,  on  a  super- 

L  D 

cavitating  hydrofoil  will  be  obtained  in  the  following 
manner.  Defining  and  CD  as  follows, 


chord 

(P  -  P, 


'/ 


dx 


chord 


(P  “  P  ) 

C-  =  /  - f-  dy 

Jjcq 


Or  we  can  write 
b 


1 
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where 


1  +  j 

P  -  P 


c  _ 


C  *5pq. 


1-  _3_  =  i_ 

1  2  l 

% 


exp  jlm(ii)  (C))| 


Im 


(C)|  =g(%)  in  Equation  (3.7), 


d x 
ds 


cosS  , 


ds  _  ,,,  e 

dT  =  -sg(-J 


-g  (£) 


y/1  +  o 


Similarly , 
b 


=  [  c  .  M-  ds1 

/  p  ds  dC  ^ 


where 

^  =  sinS . 
ds 

3 . 2  PARTIALLY  CAVITATING  FLOW  THEORY 

3.2.1  Double  Wake  Model 


One  of  the  most  difficult  problems  in  formulating  a  cavitating 
flow  theory  lies  in  construction  of  the  flow  configuration. 

The  degree  of  difficulty  increases  in  the  partially  cavitating 
fl'-'w,  particularly  for  the  flow  around  the  foil  having  a  blunt 
trailing  edge.  The  turbulent  wake  behind  the  foil  must  be 
treated  in  a  somewhat  different  manner  from  the  cavity  itself. 
In  the  present  theory  we  applied  a  "double  wake"  model  as  is 
shown  in  Figure  3.4.  In  this  model  the  wake  consists  of  two 
parts,  i.e.,  the  near  wake  and  the  far  wake.  Inside  the  far 
wake  the  dynamic  pressure  is  gradually  recovered  towards  the 
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infinity  whereas  that  of  the  near  wake  is  somewhat  in  between 
and  p  .  In  terms  of  the  velocity,  the  near  wake  wall 
velocity  is  thus  expressed  as  follows 

%  “  %  (if) +  qi{1~(ir)j  •  <3-13) 

As  was  described  in  the  report  of  Furuya  and  Maekawa  (1979), 
we  have  another  reason  for  using  the  formula  (3.13) .  When  one 
calculates  the  force  coefficients  of  cavitating  hydrofoils 
having  blunt  trailing  edges,  it  is  preferable  to  have  a  smooth 
transition  between  the  partially  cavitating  and  supercavitating 
conditions.  Equation  (3.13)  allows  us  to  smoothly  blend  the 
partially  cavitating  flow  into  the  S/C  flow.  The  near  wake 
wall  velocity  is  gradually  equated  to  qc  as  l  approaches  the 
chord  length. 

The  physical  flow  field  shown  in  Figure  3.4  is  mapped  onto  the 
potential  plane  shown  in  Figure  3.5  As  is  shown  in  this  poten¬ 
tial  plane  of  Figure  3.5,  the  velocity  q  and  flow  angle  9  on 
the  upper  and  lower  far  wake  boundaries,  are  assumed  to  have 
the  identical  values  for  the  same  potential  i> .  This  also 
requires  the  condition  that  5  at  the  upper  cavity  end  point  is 
equal  to  5  at  the  lower  cavity  end  point.  It  is  realized  that 
this  is  identical  to  the  single  spiral  vortex  model  of  Tulin 
(1964)  used  for  the  present  double  wake  flow  configuration. 

The  potential  plane  is  once  more  transferred  onto  the  upper 
half  of  the  ;-plane  with  a  mapping  function 


or 

W 


(3.14) 


20 


s 


where  denotes  the  potential  at  the  point  corresponding  to 
the  cavity  end  point  and  c  =  ia  corresponds  to  the  infinity 
in  the  physical  plane. 

The  boundary  conditions  on  the  real-axis  ~  are  now  expressed 


The  solution  for  u  will  be  obtained 
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By  taking  a  as  one  of  the  unknown  parameters,  we  have  a  total 
of  seven  unknown  quantities,  i.e.,  a,  b,  c,  f,  a,  P,  i  ,  re¬ 
quiring  seven  equations  to  determine  them  uniquely. 


Boundary  conditions  which  have  not  been  used  yet  will  be 
applied  as  follows: 

(i)  at  infinity  -  ia) 

jj(  ia)  =  i  In  A— :  2  equations 

(ii)  length  of  the  first  arc  = 

=  s^(-l)  where  in  Eqn.  (3.18)  :  1  equation 

(iii)  length  of  the  second  arc = 

S2  =  s2 where  s^  in  Eqn.  (3.22)  :  1  equation 


(iv)  cavity  end  lands  on  the  upper 
surface  of  the  blade 


1  equation 


(v)  cavity  length  matches  the  specified  length:  1  equation 


(vi)  foil-cavity-wake  system  closure  condition 


1  equation 


Application  of  these  conditions  to  Equation  (3.15)  yields  the 
following  seven  equations: 
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first  arc  =  S^ 


On  the  wetted  part  of  the  streamline, 


dz  =  ei3l 
ds^  ’ 


ds  =  e'i3l  -S-1^  —  d- 

i  qw  <*c 


For  -l<£<b,  u(c,)  can  be  written  as  follows 
ig1  U)  +  8X  U)  ,  -1  <  4  <  0 

w(5)  = 

|ig1  U)  +  8X  U)  +  *  ,  0  <  *  <b 


where 


I 


_ 1 _ 

Hz '  + 1 )  ( b  ZHEUZH 

V  £ '  -  c 


Therefore, 


( ;) 


where 


sg  ( Z ) 


£  > 


0 


<  0 


(3.18) 


and 


dW 

3T  - 


2C 

(a2  +  ?2)2  , 


(3.19) 


Thus 


f 


S.  -s. (-1)  =  0 


(3.20) 


where 


1 


» 


(vi)  foil-cavity-wake  system  closure  condition; 

For  zero  source  strength  within  the  contour  enclosing 
the  system,  we  must  have 


Re 


tj  ( ; )  dW 


=  0 


Similar  to  the  S/C  theory,  the  condition  reduces  to 

-bn.  J- 


or 


f7  =  Rejuj'  (ia)|  =  Re|HH'(X~i  u>(ia)  +  H(ia)  G'  (ia)j  =  0, 


where 

H(;) 

and 


=  JR  +  1)(S  - 


-  b)(c  ~  f) 

c 


G(^)  =  The  remaining  part  of  the  right-hand  side  of 
Equation  (3.15). 


The  lift  and  drag  coefficients  for  the  present  case  will  be 
obtained  in  a  slightly  different  way  from  that  for  the  S/C 
foil  theory.  The  lift  L  is  defined  as  follows, 


chord 

l  =y  dx 
0 

where  Pj  and  p^  denote  the  pressures  on  the  lower  and  upper 
portion  of  the  foil,  respectively.  Then,  the  normalization 
will  provide  the  lift  coefficient  CL 


1 


.25) 
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CL  " 


*sPq^  'Chord 


=  /  — ^ — dx  -/-A^ 

J  i5Pq1  -chord  J  Hoq1  -chc 


-  dx 

chord 
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Similarly,  the  drag  coefficient  CD  will  be  given 
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Like  any  other  nonlinear  theory,  it  is  a  hard  task  to  obtain 
the  first  numerically  converged  solution  for  a  posed  problem. 

Once  it  is  obtained,  the  change  of  incidence  angle,  cavity 
length  and/or  blade  profile  shape  will  not  provide  any  difficulty 
in  finding  the  solution  with  the  already  obtained  solution  used 
as  a  starting  point  for  the  iteration  loop.  The  present  partially 
cavitating  theory  with  the  double  wake  model  was  the  most  dif¬ 
ficult  theory  for  finding  the  first  converged  solution.  It 
took  us  many  months  of  computer  program  debuggings  and  trials. 

It  was  finally  found  that  Equation  (3.25)  has  a  behavior  similar 
to  (log  a)/a  with  respect  to  "a"  for  all  other  parameters  fixed. 

It  means  that  this  function  has  two  roots  for  "a",  i.e.,  1  and  *>. 
Therefore,  if  the  initial  guess  for  "a"  is  large,  Newton’s 
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method  calculates  the  slope  for  convergence  in  totally  the 
opposite  direction  and  finds  M  to  be  the  solution.  This  fact 
was  discovered  by  plotting  various  multidimensional  functions 
for  each  parameter.  By  starting  a  rather  small  value  for 
"a",  the  convergence  for  the  iteration  loop  was  immediately 
obtained.  The  computer  cost  of  calculating  the  solution  for 
one  problem  was  about  $15  with  the  7600  CDC  computer  at  the 
Lawrence  Berkeley  Laboratory. 

3.2.2  Open  Wake  Model 

This  "open  wake  model"  is  a  rather  simple  model,  as  is 
shown  in  Figure  3.7  which  is  similar  to  that  used  for  partially 
cavitating  plano-convex  cascade  flows  (see  Furuya  (1978)).  In 
this  model  it  is  assumed  that  the  pressure  behind  the  foil 
immediately  recovers  that  of  the  uniform  flow  and  continues  to 
the  downstream  infinity.  This  wake  is  considered  not  to  close 
and  thus  is  called  here  "open  wake  model". 

The  technique  of  solving  the  problem  with  this  model  is  iden¬ 
tical  to  the  previous  one,  first  mapped  onto  the  potential 
plane  W  in  Figure  3.8  and  then  C -plane  in  Figure  3.9.  Since 
this  model  and  the  previous  one  are  almost  identical  in  the 
mathematical  formulation,  only  the  differences  will  be  described 
in  the  following.  The  first  one  is  the  mapping  function  between 
the  W-plane  and  c-plane,  i.e., 

W  =  AC2,  (3.28) 

therefore  dW/dC  in  various  equations  of  the  previous  section 
should  be  replaced  by 


1=  2^. 


3.29) 


The  second  difference  is  that  the  constant  term  P  for  w  in 
Equation  (3.15)  must  be  zero.  This  condition  stems  from  the 
fact  that  w  should  be  finite  as  since  the  infinity  point 
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in  the  ';-plane  corresponds  to  the  infinity  point  in  the  physical 
plane  z.  This  point  is  quite  different  from  the  double  wake 
model  in  which  the  infinity  point  in  the  ;-plane  corresponds 
to  the  end  point  of  cavity  where  the  logarithmic  singularity 
is  allowed  to  exist. 


The  third  one  is  associated  with  the  normalization  velocity  in 

the  definition  for  w.  Since  the  wake  velocity  q  is  assumed 

2  ^w 

to  be  q^  in  the  present  model, 
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(3.30) 


where  q^  is  normalized  to  be  unity.  It  means  that  all  q^ ' s 
in  the  previous  flow  model  should  be  replaced  by  one. 


The  number  of  unknown  parameters  here  is  only  five,  i.e..  A, 
b,  c,  f,  a.  Boundary  conditions  for  determining  these  param¬ 
eters  include  (i)  through  (v)  of  the  double  wake  model  case, 
excluding  the  closure  condition.  It  must  be  mentioned  that 
only  one  condition  is  available  from  (i)  since  the  imaginary 
part  of  w  as  ;  •*  is  always  zero.  Furthermore,  a  slight  modifi¬ 
cation  was  made  here  for  the  condition  (iv)  in  Equation  (3.23). 
Instead  of  the  end  point  of  cavity  landing  right  on  the  upper 
part  of  the  blade,  we  introduced  the  turbulent  boundary  layer 
thickness  after  cavity  collapse  into  this  condition.  There¬ 
fore,  this  condition  now  reads 

f5  5  yc  -(W  *  yb  *  0 


where 


y^  =  turbulent  boundary  layer  thickness. 

This  concept  was  tried  for  the  double  wake  model.  However, 
it  had  very  little  effect  on  the  final  flow  characteristics 
determined  and  thus  was  not  explained  in  the  previous  section. 
The  selection  for  y0  is  arbitrary  and  will  only  be  determined 
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based  on  experimental  data.  This  concept  has  been  introduced 
simply  on  the  trial  basis  for  the  partially  cavitating  flow 
theory  and  will  require  a  further  refinement  for  more  rigorous 
application.  We  can  use  the  same  formula  for  the  force  coef¬ 
ficient  calculations  as  those  described  for  the  double  wake 
model  simply  by  applying  the  changes  of  Equations  (3.29)  anc 
(3.30)  . 


* 
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RESULTS 


4.1  FLOW  OBSERVATION 

Figures  4.1  through  4.5  show  various  flow  patterns  with  the 
cavitation  number  reduced  at  angles  of  incidence  0°  to  10° 
with  2  degrees  increment.  The  top  and  side  views  of  these 
photos  were  taken  simultaneously  by  using  a  strobe  light 
method.  The  duration  of  light  was  about  10  microseconds. 

The  cavitation  numbers  used  for  presenting  the  results  are 
those  corrected  according  to  the  method  mentioned  in  Section  2.4 
so  that  they  must  be  close  to  the  actual  pressure  inside  the 
cavity.  We  ensured  this  point  by  comparing  these  corrected 
cavitation  numbers  with  those  of  measurement  for  several  cases. 
It  must  be  mentioned  that  we  were  unable  to  use  the  measured 
cavity  pressures  throughout  the  present  study  since  the  data 
became  unreliable  due  to  the  air  leakage  of  the  line.  This 
experimental  problem  could  not  be  repaired  during  the  short 
period  of  tunnel  time  available  to  us. 

At  a  =  0°,  it  is  seen  from  Figure  4.1  that  the  cavity  starts 
at  the  foil  base  and  stays  there  all  the  way  until  the  tunnel 
choking  condition.  We  can  observe  several  cavitation  streaks 
of  finger-type  starting  from  the  leading  edge  at  smaller  cavi¬ 
tation  numbers  but  should  categorize  this  flow  as  a  base- 
cavitated  flow.  At  a  =  2°,  the  cavitation  now  occurred  both 
from  the  leading  edge  and  trailing  edge.  At  an  a  of  about  .35, 
these  two  cavities  merged  together  and  thus  the  flow  pattern 
became  a  supercavitating  condition.  The  f 1  >v  phenomena  observed 
at  the  time  of  merging  was  rather  abrupt.  One  of  the  original 
objectives  of  the  present  study  was  to  investigate  the  dis¬ 
crepancy  between  the  theory  and  experiment  for  a  supercavitating 
propeller.  The  discrepancy  appeared  in  a  comparison  for  the 
thrust  coefficient  at  a  higher  advance  speed  range  (see  the 
report  of  Furuya  (1978)).  This  is  the  range  of  J  at  which  the 
blade  will  experience  the  partial  cavitation  since  the  local 
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flow  incidence  angles  at  each  blade  section  become  small,  e.g., 
0°^2°.  in  the  analysis  of  Furuya  (1978)  we  assumed  a  regular, 
single  partial  cavity,  whereas  the  experiment  showed  either 
base-cavity  or  double  cavity  for  the  partial  cavity  occurring 
at  ct  =  0°~2°.  It  seems  quite  possible  that  the  discrepancy 
in  the  propeller  performance  prediction  method  is  attributable 
to  this  erroneous  perception  of  the  flow  configuration  at  the 
partially  cavitating  flow  regime. 

For  the  incidence  angle  larger  than  4°,  no  base  cavity  was 
observed  any  more  and  all  cavities  started  from  the  leading 
edge  and  grew  longer  with  decrease  in  pressure.  Most  of  the 
cavities  observed  here  were  of  bubble  type.  The  cavity  walls 
became  clear  and  glassy  for  the  cases  of  high  incidence  angle 
with  long  cavities  (see  Figures  4.4  (f)  and  4.5  (f)). 

Although  we  did  not  observe  any  natural  instability  near 
^c/chord  =  1  such  as  that  reported  in  the  paper  of  Wade  and 
Acosta  (1967) ,  we  had  difficulty  in  keeping  the  flow  steady 
at  a  desired  cavity  length  for  large  incidence  angles  during 
the  partially  cavitating  condition.  The  phenomena  will  clearly 
be  seen  in  Figures  4.5(c)  and  (d)  ;  immediately  after  7c/chord  = 
0.5  was  achieved,  the  cavity  extended  abruptly  beyond  the 
chord  length  with  a  slight  decrease  in  pressure  and  could  not 
hold  at  any  point  in  between  them.  It  is  for  this  reason  that 
only  a  few  data  points  were  obtained  at  the  p/c  condition  for 
a  =  8°  and  10°,  as  will  be  seen  in  the  force  data. 

4.2  FORCE  COEFFICIENTS 

Data  corrections  for  the  measured  lift  and  drag  coefficients 
were  made  in  the  way  explained  in  Section  2.4  and  they  are 
shown  in  Figures  4.6  and  4.7  at  .n  =  0°  through  10°.  The 
cavitation  numbers  used  here  are  again  the  ones  with  correc¬ 
tions  and  thus  are  supposed  to  be  a  .  .  's. 

cavity 
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All  data  points  taken  during  the  experiment  are  shown  in  these 
figures  and  not  a  single  point  is  excluded.  Figures  4.8  and 
4.9  show  composite  lift  and  drag  coefficients  made  by  inter¬ 
polating  the  experimental  data.  It  is  seen  from  these  figures 
that  all  lift  curves  show  a  similar  trend  as  a  function  of  J 
except  for  a  =  0°  as  mentioned  above.  It  must  be  noted,  how¬ 
ever,  that  the  hump  near  the  transition  point,  i.e.,  between 
the  S/C  and  P/C  regimes,  is  the  largest  at  a  =  2°  and  gradually 
becomes  unnoticeable  as  a  increases. 

For  the  supercavitating  range,  particularly  at  small  cavitation 
numbers,  the  lift  forces  increase  with  the  incidence  angle  but 
this  relationship  overturns  for  short  cavity  cases.  This  is 
the  behavior  which  all  theories,  including  nonlinear  theories, 
fail  to  predict  since  the  effect  of  cavity  collapse  and  sub¬ 
sequent  turbulent  mixing  near  the  foil  cannot  be  accurately 
represented  mathematically. 

The  newly  developed  theories,  one  for  the  supercavitating  flows 
and  two  for  the  partially  cavitating  flows,  were  used  for  com¬ 
parison  with  the  experimental  data.  The  results  are  shown 
only  for  the  incidence  angles  of  4 ,  6  and  8  degrees  in  Figures 
4.6  and  4.7.  The  supercavitating  flow  theory  for  the  lift 
coefficient  agrees  well  with  the  experiment  for  small  cavita¬ 
tion  numbers  at  which  the  cavity  lengths  are  long,  but  under¬ 
estimates  them  for  short  cavity  cases  .  For  the  drag  coefficients 
the  overall  prediction  capability  of  the  theory  is  reasonable, 
however,  no  definite  trend  for  the  discrepancy  exists.  This 
fact  seems  to  confirm  the  difficulty  in  accurate  measurements 
and  data  reduction  for  the  drag  force  in  this  type  of  experiment. 

For  the  partially  cavitating  regime,  the  two  new  theories  pro¬ 
vide  analytical  data  for  the  lift  and  drag  forces.  The  double 
wake  model  substantially  underpredicted  both  forces,  whereas 
the  open  wake  model  showed  an  accurate  trend  with  a  choice 
of  the  wake  thickness  yg  in  Equation  (3.31)  as  shown  in  Figure 


35 


I 


4.10.  As  has  been  mentioned  earlier,  it  remains  to  be  seen 
whether  or  not  this  choice  of  in  Figure  4.10  can  be  applied 
to  the  other  types  of  foils  for  accurate  prediction  of  lift 
forces.  It  must  be  mentioned,  however,  that  the  theory  showed 
a  numerical  instability  as  £  /c  becomes  greater  than  0.8, 
similar  to  that  reported  in  Furuya  (1980) . 

It  was  a  disappointment  that  the  double  wake  model  failed  to 
predict  the  forces  accurately.  As  will  be  seen  later,  the 
only  problem  of  this  method  stems  from  the  poor  capability  for 
predicting  the  cavitation  number  against  the  specified  cavity 
length.  If  this  were  done  properly,  the  force  coefficients 
would  be  predicted  accurately.  It  can  readily  be  seen  from  the 
fact  that  the  force  coefficient  for  the  partial  cavity  flow 
will  be  determined  mainly  by  the  cavity  pressure.  Efforts  were 
made  on  improving  the  theory.  These  included  the  wake  pressure 
change  by  using  Uc/c)n  in  Equation  (3.13)  with  n  =  2  and  3, 
as  well  as  the  application  of  the  wake  thickness  y^ ,  similar 
to  that  used  in  the  open  wake  model,  but  without  success. 

It  seems  that  more  investigation  will  be  worthwhile  since  the 
method  is  physically  more  reasonable  and  predicted  the  drag 
forces  better  than  the  infinite  wake  model.  Finally,  it  must 
be  mentioned  for  the  limiting  case  as  i  *  0,  the  two  theories 
matched  well. 

4 . 3  CAVITY  LENGTH  vs .  CAVITATION  NUMBER 

The  cavity  length  as  a  function  of  cavitation  number  is  one 
of  the  important  features  both  in  design  and  analysis  work  for 
high  speed  foils  and  propellers.  Figures  4.11(a)  through  (e) 
show  such  data  in  comparison  with  some  theoretical  data. 

Except  for  the  case  of  a  =  4°,  the  supercavitating  foil  theory 
predicts  this  relationship  poorly.  Both  the  experiment  and 
theory  are  to  blame  for  the  discrepancy.  For  the  short  cavity 
cases,  the  theory  seems  to  fail  to  represent  the  real  flow 
pattern  accurately,  whereas  for  the  long  cavity  cases,  the 
water  tunnel  chokes  so  that  the  tunnel  wall  effect  comes  into 
the  picture  and  causes  an  error  in  measurement. 
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As  for  the  partially  cavitating  regime,  the  results  made  with 
the  open  wake  model  match  well  with  the  experimental  data, 
whereas  those  of  the  double  wake  model  have  a  substantial  dis¬ 
crepancy.  It  must  be  mentioned  again  that  success  of  any 
partially  cavitating  flow  theory  will  require  an  accurate 
prediction  of  the  cavitation  number  for  the  given  cavity 
length. 
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CONCLUSION 


As  a  preparatory  phase  for  the  cascade  experiment  to  be  con¬ 
ducted  in  the  following  phase,  single-foil  experiments  were 
carried  out.  The  blade  profile  shape  of  the  foil  selected  here 
is  that  of  Hydronautics  7607.02  supercavitating  propeller  at 
50%  radial  station.  It  consisted  of  a  modified  Tulin  two-term 
camber  with  a  blunt  trailing  edge,  which  is  considered  to  be 
a  typical  supercavitating  foil  profile  shape.  The  foil  was 
installed  in  the  two-dimensional  test  section  of  the  High  Speed 
Water  Tunnel  at  the  California  Institute  of  Technology  and 
tested  over  a  full  range  of  cavitation  numbers.  Since  it  was 
found  that  there  existed  no  appropriate  nonlinear  theories 
for  comparison,  the  major  effort  after  the  experiment  was 
directed  towards  the  development  of  such  theories.  We  have 
developed  three  nonlinear  cavitv-flow  theories,  one  for 
supercavitating  flows  and  two  for  partially  cavitating  flows. 

It  is  believed  that  these  new  nonlinear  theories  will 
provide  tools  for  design  and  analysis  of  cavitating  foils 
and  propellers. 

The  success  of  the  present  experiments  for  the  single  foil  has 
led  us  to  design  dummy  foils  which  will  form  a  cascade  con¬ 
figuration  for  the  following  cascade  experiments.  Four  dummy 
blades  will  be  constructed  for  achieving  the  solidity  of  unity, 
but  two  dummy  blades  with  the  center  blade  may  be  tested  for 
the  solidity  of  0.5.  The  drawing  for  the  dummy  blades  has 
already  been  made  for  fabrication.  The  algorithm  for  HP9845 
of  controlling  the  downstream  wall  angle  for  the  cascade 
experiment  was  also  developed  during  the  present  work. 
Incidentally,  this  HP9845  was  used  even  for  the  single  foil 
experiments  as  an  on-site  data  reduction  system  to  help  us 
reduce  the  experimental  errors  or  discover  any  abnormalities 
during  the  work. 

The  conclusions  drawn  from  the  present  preliminary  work  are 
summarized  as  follows: 
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1)  The  foil  used  for  the  present  experiment  could  not 
maintain  the  leading  edge  cavitation  at  low  incidence 
angles  such  as  2°  or  less.  At  2°  the  cavity  appeared 
both  from  the  leading  edge  and  trailing  edge  and  by 
reducing  it  to  0°,  the  cavity  stays  at  the  trailing 
edge  over  a  complete  range  of  cavitation  number. 

2)  Due  to  the  flow  pattern  change  from  large  to  small 
incidence  angles,  the  lift  and  drag  force  curve  trends 
also  changed.  It  seems  that  this  was  not  predicted 
during  the  design  work  and  it  may  be  for  this  reason 
that  the  thrust  coefficient  of  7607.02  S/C  propeller 
was  overpredicted  (see  the  paper  of  Peck  (1977)  for 

the  comparison  between  the  design  value  and  experimental 
data).  This  propeller  was  designed  under  the  assumption 
that  the  leading  edge  cavity  would  exist  at  local  flow 
incidence  angles  of  2°,  but  the  assumption  seems  incor¬ 
rect. 

3)  For  exactly  the  same  reason  as  above,  the  overprediction 
of  the  thrust  at  the  large  J  range  with  the  off-design 
prediction  method  of  Furuya  (1978)  may  have  arisen.  We 
incorrectly  used  the  force  data  which  belonged  to  the 
leading  edge  cavity  flows  even  for  =  2°.  For  such 
low  incidence  angles,  we  should  have  applied  the  data 
for  base  cavity  flows.  The  forces  there  are  much 
smaller  than  those  of  the  S/C  condition,  as  has  been 
seen  in  the  experimental  results. 

4)  It  is  an  extremely  difficult  task  to  determine  theo¬ 
retically  at  which  incidence  angle  the  cavity  flow 
pattern  changes  from  the  leading  edge  cavity  to  base 
cavity,  particularly  for  supercavitating  foils  having 
blunt  trailing  edges.  This  is  one  of  the  most  impor¬ 
tant  features  for  the  design  of  supercavitating  foils 
and  propellers  to  determine  their  performance  accurately 
and  thus  will  require  further  investigation. 
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5)  The  supercavitating  flow  theory  with  the  single  spiral 
vortex  model  can  accurately  predict  the  force  coef¬ 
ficients  for  the  cases  of  cavity  length  larger  than 
1.5  chords,  but  underpredicts  for  shorter  cavity  cases. 

6)  The  partially  cavitating  flow  theory  with  the  double 
wake  model  failed  to  predict  the  right  cavitation 
number  for  the  specified  cavity  length,  thus  substan¬ 
tially  underpredicted  the  force  coefficients. 

7)  The  partially  cavitating  flow  theory  with  the  open 
wake  model  could  fit  the  lift  curves  with  the  experi¬ 
mental  data  by  adjusting  the  turbulent  mixing  wake 
thickness.  The  wake  thickness  was  determined  as  a 
function  of  a  for  the  present  cases  but  needs  justi¬ 
fication  for  its  applicability  to  other  foil  profiles. 

8)  Both  partially  cavitating  flow  theories  predicted  the 
drag  forces  with  poor  accuracy. 

9)  The  prediction  accuracy  of  all  the  theories  developed 
here  degrades  as  soon  as  the  cavity  collapses  near  the 
trailing  edge  of  the  foil  both  for  the  supercavitating 
and  partial  cavitating  conditions.  It  seems  that  highly 
turbulent  flow  phenomena  near  the  cavity  collapse  will 
need  a  better  flow  modeling  for  mathematical  formulation. 

10)  The  analytical  prediction  capability  for  the  partial 

cavity  length  is  important  for  the  high  speed  propeller 
analysis.  The  theory  with  the  infinite  wake  model 
showed  a  promising  feature  for  this  purpose  but  its 
universal  application  still  remains  to  be  seen. 
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TABLE  2.1 


MODEL  BLADE  OFF-SET  POINTS 

(Hydronautics  7607.02  S/C  Propeller  at  50%  Radial  Station) 


STATION 

OFFSET 

INCHES 

TOTAL 

THICKNESS 

Percent 

Chord 

Distance 
from  L.E. 

Upper 

(Back) 

Lower 

Face 

Inches 

Percent 

0 

0 

0 

0 

0 

0 

5 

.160 

.023 

-.009 

.032 

1.0 

10 

.320 

.039 

-.016 

.055 

1.7 

15 

.480 

.051 

-.023 

.074 

2.3 

20 

.640 

.060 

-.029 

.089 

2.8 

25 

.800 

.063 

-.036 

.104 

3.3 

30 

.960 

.075 

-.044 

.119 

3.7 

35 

1.120 

.082 

-.052 

.134 

4.2 

40 

1.280 

.088 

-.060 

.148 

4.6 

45 

1.440 

.093 

-.068 

.161 

5.0 

50 

1.600 

.098 

-.077 

.175 

5.5 

55 

1.760 

.10  3 

-.086 

.189 

5.9 

60 

1.920 

.107 

-.096 

.203 

6.3 

65 

2.080 

.110 

-.107 

.217 

6.8 

70 

2.240 

.114 

-.117 

.231 

7.2 

75 

2.400 

.116 

-.129 

.245 

7.7 

80 

2.560 

.119 

-.140 

.259 

8.1 

85 

2.720 

.120 

-.153 

.273 

8.5 

90 

2.880 

.122 

-.167 

.289 

9.0 

95 

3.040 

.122 

-.183 

.305 

9.5 

100 

3.200 

.121 

.323 

10.1 

FIGURE  3.1 

FLOW  CONFIGURATION  FOR  S/C  FOIL 


FIGURE  3.2 
POTENTIAL  PIANE 
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FIGURE  3.3 
TRANSFORM  PLANE 


FIGURE  3.4 

FLOW  CONFIGURATION  OF  DOUBLE  WAKE 
MODEL  FOR  PARTIALLY  CAVITATING  FOIL 
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FIGURE  3.5 
POTENTIAL  PLANE 


FIGURE  3.6 
TRANSFORM  PLANE 
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FOR  PARTIALLY  CAVITATING  FOIL 


FIGURE 


FIGURE  4.2  SIMULTANEOUS  VIEWS  FROM  TOP  AND  SIDE  FOR  INCIDENCE  ANGLE 


FIGURE  4.3  CONTINUED 


NEOUS  VIEWS  FROM  TOP  AND  SIDE  FOR  INCIDENCE  ANGLE 


00 

(N 

II 

D 


V  ■'*, 


:  *-H%& 


i 


5  3 


FIGURE  4.4  CONTINUED 
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FIGURE  4.5  CONTINUED 
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FIGURE  4.6  CONTINUED 


CONTINUED 


FIGURE  4.7  CONTINUED 
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1NUED 


CAVITATION  NUMBER 


Wake  Thickness  y  /chord 
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Incidence  Angle  a  (degrees) 
FIGURE  4.10 

WAKE  THICKNESS  IN  EQUATION  (3.31) 
USED  FOR  THE  INFINITE  WAKE  MODEL 
AT  lc/c  =0.5  WITH  A 
PROPORTIONALITY  LAW  FOR  CHANGE  OF  1 
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CAVITATK  .  NUMBER 
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FIGURE  4.11  CONTINUED 


